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Abstract 

A novel flow state consisting of two oppositely travelling waves (TWs) with oscillating amplitudes 
has been found in the counterrotating Taylor-Couette system by full numerical simulations. This 
structure bifurcates out of axially standing waves that are nonlinear superpositions of left and right 
handed spiral vortex waves with equal time-independent amplitudes. Beyond a critical driving the 
two spiral TW modes start to oscillate in counterphase due to a Hopf bifurcation. The trigger for 
this bifurcation is provided by a nonlinearly excited mode of different symmetry than the spiral 
TWs. A three-mode coupled amplitude equation model is presented that captures this bifurcation 
scenario. The mode-coupling between two symmetry degenerate critical modes and a nonlinearly 
excited one that is contained in the model can be expected to occur in other structure forming 
systems as well. 
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Many nonlinear structure forming systems that are driven out of equilibrium show a 
transition to traveling waves (TWs) as a result of an oscillatory instability 1]. In the 
presence of spatial inversion symmetry in one or more directions also a standing wave (SW) 
solution bifurcates which is a nonlinear superposition of the two symmetry degenerated, 
oppositely propagating TWs with equal amplitudes. 

SWs and TWs have a common onset as a result of a primary Hopf bifurcation. But at 
onset only one of them is stable 2|, y|. Furthermore, there are mixed patterns with non-equal 
amplitude combinations of the degenerated TWs that arise, e.g., via secondary bifurcations 
at larger driving. The variety with temporally constant, non-equal TW amplitudes can 
provide a stability transferring connection between TWs that, e.g., are stable at onset and 
SWs that become stable later on 4j. 

The variety in which the TW amplitudes oscillate in time is the subject of this paper. 
This solution bifurcates out of the SW via a Hopf bifurcation. To be concrete we investigate 
wave structures consisting of spiral vortices in the annular gap between counter rotating 

n n 

concentric cylinders of the Taylor- Couette system [5|, |GJ. To that end we have performed 
numerical simulations of the Navier-Stokes equations (NSE) to reveal the bifurcation prop- 
erties as well as the spatiotemporal structure of the novel oscillating mixed wave states. In 
addition, we provide coupled three-mode amplitude equations that capture this bifurcation 
to explain the underlying mode-coupling mechanism. We are not aware that these states 
have been reported so far in the Taylor- Couette literature. Furthermore, one can expect 
that the mode-coupling mechanism between two symmetry degenerate critical modes and 
the nonlinearly excited one that is described by our coupled amplitude equations and that 
drives the oscillatory instability is operating in other pattern forming systems as well. 

The waves are realized by left handed spiral vortex (L-SPI) and right handed spiral vortex 
(R-SPI) structures that are mirror images of each other. The azimuthal advection by the 
basic circular Couette flow (CCF) rotates both like rigid objects into the same direction 
as the inner cylinder [?]]. As a result of the enforced rotation the phases of L-SPI and R- 
SPI travel axially into opposite directions. This system offers an easy experimental and 
numerical access to forward bifurcating TWs and SWs that are called ribbons (RIBs) 8| in 
the Taylor-Couette literature. Being a nonlinear superposition of L-SPI and R-SPI the RIB 
structure also rotates azimuthally, however, such that its oscillations in axial direction form 
a SW. 
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Here we elucidate how such stable SWs loose stability to an oscillating state via a Hopf 
bifurcation. Therein, the interaction with another nonlinearly excited, non-traveling mode 
induces the TW constituents of the SW, i.e., the L-SPI and the R-SPI to oscillate in counter- 
phase around a common mean. These oscillating mixed wave states that we call oscillating 
cross spirals (O-CR-SPI) are quite robust. Thus, they should easily be observable in exper- 
iments. 

All these spiral structures are axially and azimuthally periodic. We have focussed our 
simulations on patterns with axial wavelength A = 1.3 measured in units of the gapwidth 
and azimuthal wave number M = 2. The numerical solutions of the NSE were obtained for 
a system with radius ratio rj = 1/2 by methods described in Q. 

Control- and order parameters - The rotational velocities of the inner and outer cylinders 
are measured by the respective Reynolds numbers R\ and R 2 . We fix Ri = 240 [12] and we 
introduce the reduced distance /x = (R% — R^)/\R^\ from the common onset of SPI and RIB 
flow at i?2 = —605.5 as control parameter. We characterize the spatiotemporal properties 
of the vortex waves using the Fourier decomposition 

f(r, tp, z,t) = J2 fm,n(r, t) e^ +nk ^ (1) 



in azimuthal and axial direction. Here one has f- m - n = f m , n with the overbar denoting 
complex conjugation. Order parameters are the moduli \A\, \B\, \C\ and the time derivatives 
Qai@b-,Qc °f the phases of the dominant modes in the decomposition (Tj[|) of, say, the radial 
velocity u at midgap u 2 ,i = A = \A\e~ teA , u 2 -i = B = \B\e~ ldB , and u 0:2 = C = \C\e~ l6c . 
Here, A and B are the amplitudes of the marginal L- and R-SPI modes. When both are finite 
as, e.g., in the SW of the RIB state their nonlinear coupling generates the m = C-mode 
below its threshold for linear growth: pure m = stationary Taylor vortices bifurcate out 
of the CCF only later on. Although |C| itself remains small compared to \A\, \B\ in the RIB 
state its feedback on A, B triggers the Hopf bifurcation of the O-CR-SPI: the oscillations 
of, say, \A\ are driven by bilinear mode couplings of BC as indicated in Fig. [TJ 
We also use the combined order parameters 

I 4| 2 i I D|2 I 4 12 _ I D|2 

s = |g| +|g| D = |g| ^ = 9c + 9B _ eA _ 7l (2) 

that are better suited to describe the bifurcation of the O-CR-SPI with oscillating D(t) and 
$(t) out of the RIB state, D = = 
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Bifurcation sequence - The pure TW shown by circles in Fig. [2] and the SW solution 
(A = B, C 7^ 0) marked by diamonds bifurcate at fi = out of the unstructured CCF. 
Initially, the SPI is stable and the RIB is unstable. But then there appears a stable cross- 
spiral (CR-SPI) solution [triangles in Fig. [2(b)] which transfers stability from the SPI to the 
RIB. The moduli and phase velocities of these three structures are time- independent. At \ih 
in Fig. [2] the RIB lose stability in a supercritical Hopf bifurcation to the novel modulated 
state of O-CR-SPI. Increasing fi further beyond the range shown in Fig. [2] the O-CR-SPI 
loses stability at R2 ~ —543 to oscillating structures with azimuthal wave number M = 1 
that are not discussed here. 

Dynamics of the modulated SW - Figure [3] shows the temporal variation of characteristic 
quantities of the O-CR-SPI over one modulation period r. Thick lines refer to /1 immediately 
above onset Thin lines show behavior at a larger value fi > (arrow in Fig. [2]) that is close 
to the end of the existence interval of O-CR-SPI. The moduli \A\, \B\ in Fig. [3(a) and 
the phase velocities 9a, Ob in Fig. [3(c) oscillate each in counter phase around a respective 
common mean. Also 6c oscillates. Furthermore, \C\ and the combined order parameter S 
show small amplitude oscillations with twice the frequency of the other quantities. Close to 
onset all oscillations are harmonic with \C\ and S being practically constant. But at yU> the 
oscillations of \A\, \B\ and 6 a, 63,^ have become quite anharmonic, whereas those of S, D, 
\C\, and 6c are still harmonic. The Fourier spectra in Figs. 3(a)- (g) of the temporal profiles 
shown by thin lines in Figs. [3(a)- (g) reflect this behavior at /i>. 

In the RIB state the phases are such that 6c it) + 0b(£) — = tt- But as a consequence 
of the Hopf bifurcation D as well as $ oscillate in the O-CR-SPI. The squares of their 
oscillation amplitudes, D 2 and $ 2 , increase at onset linearly with /i with a subsequent 
quadratic correction, cf. Figs. [5(c)- (d). The monotonous decrease of the modulation period 
r is shown in Fig. [5(a). Note that the modulation amplitudes of S in Fig. [5(b) and also of 
\C\ remain very small compared to those of D and $. 

Amplitude equations - The Hopf bifurcation behavior and the dynamics close to the 
transition from RIB to O-CR-SPI can be explained and described within a three-mode 
amplitude-equation approach. It reveals (i) how the rotationally symmetric C-mode is 
generated nonlinearly via the interaction of A and B, i.e., of the M=2 SPI constituents in 
the RIB and (ii) how then C - after it has reached a critical size beyond fin ~ induces 
amplitude oscillations in A and B. 
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Invariance under axial translation and reflection of the Taylor-Couette system [d] restricts 
the form of the three coupled amplitude equations to 

A = AG (\A\ 2 , \B\ 2 , \C\ 2 ) + ikBC, (3a) 
B = B G (\A\ 2 , \B\ 2 , \C\ 2 ) + ikAC, (3b) 
C = C H (|A| 2 , |5| 2 , |Cf ) + k AB. (3c) 

With G{\A\ 2 ,\B\ 2 ,\C\ 2 ) = G(\B\ 2 ,\A\ 2 ,\C\ 2 ) and H(\A\ 2 , \B\ 2 , \C\ 2 ) = H(\B\ 2 ,\A\ 2 ,\C\ 2 ) 
the equations are invariant under the operation (A, B, C) (B, A, C) which reflects the 
axial inversion symmetry z <-> —z. The functions G = G + %G" and H = H + iH" are 
complex. The superscripts ' and " identify the real and imaginary parts, respectively. The 
coupling constants k and k,q are real. 

Since only invariance under translation and reflection along one spatial direction has been 
used in deriving Eqs. ([3]) our description of the phenomenon of a SW with oscillating TW 
components in terms of Eqs. (jHJ) potentially applies to all bifurcating systems with 0(2) 
symmetry in the center manifold, which is quite common. 

In the following we discard the coupling term n AB. It is small in our case and, more 
importantly, we checked that it is not relevant for driving the Hopf oscillations. They are 
generated by the coupling terms in ( |3a|) and ( I3bl) as we shall show in the next section. 

The mechanism causing the Hopf bifurcation into the modulated SW can be better iso- 
lated by rewriting the amplitude equations ([3]) in terms of the combined order parameters 

3 = 2[D G'_ + S G' + ], \C\ = |C| H', (4a) 
D = 2 [S G_ + D G' + ] - 2k\C\S* sin $, (4b) 
6 = 2G"_ + 2k^Dcos§-H'\ (4c) 

where S* = Sy/1 - (D/S) 2 ~ S. Here we defined G± = (G ± G)/2. Note that G + and H' 
(G- and H ) are even (odd) in D result of the inversion symmetry. Hence, eqs. 

(l4aj) are even in D. This in turn explains that S and \C\ oscillate with twice the frequency 
of the other quantities in Fig. On the other hand, Eq. (I4"cl) is odd in D and causes the 
absence of a peak at 2/r in Fig. IU(g). 

We have determined the specific functions G and H for our specific system via fits to the 
numerically obtained bifurcation branches of SPI, RIB, and CR-SPI in Fig. [2] and to the 
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pure Taylor vortex solution (not shown here) with A = = B, C ^ 0, and half the spiral 
wavelength. This produces the bifurcation behavior close to the Hopf threshold well. Note, 
however, that the Hopf bifurcation is a universal phenomena of systems like (131^0) that is 
not specific to the Taylor- Couette system. This is most easily understood with the help of 
the universal small-/} expansion of (jlj) that results from the symmetry properties. 

Hopf bifurcation - For small D, i.e., close to the Hopf bifurcation threshold we can use 
the expansions 

G + = Gf + O (D 2 ) , G_ = G W D + O (D 3 ) , (5a) 
H' = + O (D 2 ) , H" = H"<»D + O (D 3 ) . (5b) 

Here the leading order terms G+\ G^\ H'(°\ H"^ still depend on S and \C\ 2 . Inserting (jSJ) 
into (j3]) and using the smallness of $ yields a simplified model that is linear in D 

S = 2S Gf\ \C\ = |C| H' {0 \ (6a) 



D 



S G w + G'i 0) 



2g"^ + 2k M _ H "W 



D-2k\C\S$, (6b) 
s - ,D. (6c) 

It explains the Hopf bifurcation out of the RIB state and the O-CR-SPI properties close to 
onset. For example, S and |C| are virtually constant because they are decoupled from D 
and $ in the model eqs. Furthermore, Eq. (16"c|) shows that $ is enslaved by D and that 
the phase shift between them is r/4 as to be seen in Fig. [3] close to This justifies the 
solution ansatz 

D(t) = Dcos(u H t), $(t) = ®sm(uj H t) (7) 

where ujh is the Hopf frequency. 

The latter is identified together with the bifurcation threshold \xu by a linear stability 
analysis of the RIB fixed point D = = $, S = SjubW, C = Crib{h) for which G+ = 
according to Eq. (|6al) . Thus, the linearized equations for the stability- relevant deviations 
from this fixed point read D = aD + 6$, $ = cD, with coefficients a = 2SG1 , b — 
—2k\C\S, c = 2G}^ + 2k-^- — H"W to be taken at the RIB fixed point. Consequently, the 
location of the zero in a(/i) determines \xu and the imaginary part of the eigenvalue at fin, 
i.e., the Hopf frequency is then given by uj 2 h = —be oc k + h.o.t, revealing that the coupling 
terms in (l3aj) and (l3bl) cause the Hopf bifurcation. Furthermore, a(/z) = a(/i — with 
positive a to ensure decay of oscillations below \i H and growth above it. 
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Conclusion - The bifurcation of a novel spiral vortex structure with oscillating TW am- 
plitudes out of an SW is shown to be triggered by the coupling to a nonlinearly excited mode 
when the latter exceeds a critical strength. Since this novel state, in which the TW ampli- 
tudes oscillate in counterphase around a common mean occurs quite robustly in a relatively 
wide parameter range it should be easily accessible to experiments. 

Our results have been obtained by full numerical simulations and explained and confirmed 
by a coupled amplitude equation model that captures the mode-coupling between two sym- 
metry degenerate critical modes and a nonlinearly excited one. Our bifurcation scenario can 
occur in all systems with an 0(2) symmetric center manifold, arising for example in systems 
with translation and inversion symmetry, which is a quite general one. It has therefor the 
potential to occur also in other structure forming systems, say, in hydrodynamics, chemical 
reactions, or biological systems etc. where any two symmetry degenerate basic modes A and 
B couple similarly to a third one, C, that is nonlinearly excited by them and that destroys 
the A = B state once C has reached a critical size. 

This work was supported by the Deutsche Forschungsgemeinschaft. 
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FIG. 1: (Color online) (a) Dominant modes and their complex conjugates in the Fourier space of 
Eq. ([1]). (b) Bilinear coupling of modes B and C (dashed arrows) that drive oscillations of mode 
A (solid arrow). 
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FIG. 2: (Color online) Bifurcation diagrams of SPI (red circles), RIB (blue diamonds), CR-SPI 
(purple triangles), and O-CR-SPI (mangenta lines and crosses) obtained from numerical simulations 
of the NSE versus fi and R2- SPI and CR-SPI are displayed only in (a) and (b); the latter shows the 
blow-up of the rectangle near the origin of (a). Shown are the squared mode amplitudes |A| 2 , |i?| 2 
(a), (b) and phase velocities 6 a, Ob (c) of the marginal modes and the same for the nonlinear 
excited mode C in (d) and (e). Filled (open) symbols denote stable (unstable) solutions with time- 
independent amplitudes. Crosses refer to temporal averages of the O-CR-SPI. Upper (lower) line 
shows the maximum (minimum) of the oscillation range indicated by vertical lines. The arrow at 
fJ-H {R2 = —587) marks the Hopf bifurcation of the modulated SWs. The second arrow at fj, = /i> 
(i?2 = —546) is inserted for later reference. 
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FIG. 3: (Color online) Time variation of O-CR-SPI over one modulation period r. The left column 
shows moduli and phase velocities of A, B, and C. In (a) and (c) solid lines refer to A and dashed 
ones to B. The right column contains the order parameters S, D, and <I> ([2]). Thick lines are 
modulation profiles close to the Hopf threshold fin and thin ones those at the larger ^> identified 
by second arrow in Fig. [2j 
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FIG. 4: (Color online) Fourier spectra of the modulation profiles shown by thin lines in Fig. [3] for 
/i = fiy (cf. arrow in Fig. [2]). Note that \C\ and S oscillates with twice the frequency of the other 
quantities and that the spectra of \C\, S, and D practically do not contain higher harmonics. The 
spectra of D, and 9c (|C| and S) contain only peaks at (21 + l)/r (21/t) with I = 0, 1,2 
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FIG. 5: (Color online) Bifurcation properties of RIB (blue diamonds) and O-CR-SPI (magenta 
lines and squares) obtained from numerical solutions of the NSE as functions of \i and Ri'- (a) 
oscillation period r of the modulation, say, of the moduli |^4| and \B\ of the O-CR-SPI, (b) S 
(thin lines delimit the oscillation range indicated by vertical bars), (c) and (d) squared oscillation 
amplitudes D of D and of respectively. 



